Small polarons and c-axis transport in highly anisotropic metals 
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Motivated by the anomalous c-axis transport properties of the quasi two-dimensional metal, 
Sr2Ru04, and related compounds, we have studied the interlayer hopping of single electrons that 
are coupled strongly to c-axis bosons. We find a c-axis resistivity that reflects the in-plane electronic 
scattering in the low and very high temperature limits (relative to the characteristic temperature of 
the boson Tboson)- For temperatures near the Xbosom a broad maximum in the resistivity can appear 
for sufficiently strong electron-boson coupling. This feature may account for the observed "metallic 
to non-metallic crossover" seen in these layered oxides, where the boson may be a phonon. 



I. INTRODUCTION 



Since the discovery of the cuprates, many studies have 
focused on the putative two dimensional (2d) non Fermi 
liquid state of the .copper oxide layer in the normal state 
of these materials^ However, the superconducting state 
is three dimensional and, at least near to the super- 
conducting critical temperature (T c ), some aspects of 
the physics-iof coupling between the planes must become 
important, a 

The transport properties in the weakly-coupled dimen- 
sion (denoted the c-axis here) are also interesting in their 
own rightu. In the normal state of many of these layered 
copper oxides, the resistivity ratio, p a b/Pc is of order 10 3 
to 10 5 but intriguingly, this ratio is temperature depen- 
dent: p c is non-metallic near T c , p c ~ T -7 , < 7 < 
2, while the in-plane resistivity has the marginal Fermi 
liquid behavior p a i, ~ T (at least at optimal dopingp. 

This dichotomous behavior cannot easily be accommo- 
dated within an anisotropic three dimensional (3d) Fermi 
surface picture. We begin with the assumption that the 
c-axis coupling is simply interlayer hopping of single elec- 
tron: 



H c = t, 



£4- 



h.c. 



(1) 



where n is the layer index. Then, at temperatures larger 
than the c-axis electron hopping matrix element t c , sim- 
ple perturbation theory (to lowest order in t c ) with the 
Kubo formula gives the c-axis conductivity to be pro- 
portional to the polarization bubble (Fig. |l|). In the fig- 
ure, the thick lines represent the renormalised in-plane 
Green's function. Hence the c-axis conductivity directly 
reflects the in-plane spectral weight. But if the in-plane 
state is a metal, then the c-axis conductivity should also 
be metallic. 

There are at least two ways to bypass this argument. 
Either the in-plane physics is unusual or the interlayer 
coupling is unconventional. Considering the first possi- 
bility, there might be a non Fermi liquid in-plane Green's 
functional. More conventionally, Ioffe and MillisQ ar- 
gued for superconducting phase fluctuations, and oth- 
ers have just assumed some phenomenological form of 
in-plane spectral weight which has strong in-plane scat- 
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FIG. 1: The order t c contribution to the c-axis conductivity. 
The thick lines represent the renormalised in-plane Green's 
function. 

tering built-i nl. In such schemes, the interlayer trans- 
port directly probes the anomalous in-plane scattering. 
Within the second possibility, Rojo and LevinEl consid- 
ered static disorder, and also boson-assisted hopping in 
the c-axis direction, while Turlakov and Leggettfj studied 
interplane and in-plane charge fluctuations. These ap- 
proaches can lead to incoherent c-axis transport, without 
the need to invoke strongly anomalous in-plane behavior. 

Motivated by experiments in a related system, the 
strontium ruthenate family, we follow the second route 
in this paper and study a boson-assisted hopping mech- 
anism for the unusual c-axis transport in the ruthenate 
systems. 

In contrast to previous work on boson-assisted 
hoppingB, in this report we treat the strong c-axis 
electron-boson coupling exactly via a canonical trans- 
formation, and then add the weak c-axis hopping per- 
turbatively. Xhus the essential physics is that of the 
small polaronE3'M This can be seen most easily after the 
canonical transformation: when the electron hops from 
one plane n + 1 to another at n, a cloud of bosons are 
created or destroyed: 



f canon. + 

cJjCn+i — ► c^c n+ iexp 

transf . 



^A n {q){a q -aL q ) 



L q 



(2) 



(The precise form of the amplitude A n (q) is shown in 
Eq. [7| below.) At temperatures lower than the character- 
istic energy of the boson, Tbosom electrons can hop from 
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one plane to another, but this hopping is suppressed due 
to the "dragging" of the boson cloud. With increasing 
T, more and more bosons are created and destroyed; due 
to this inelastic scattering process, the hopping electron 
acquires a large imaginary part in the self-energy: in- 
ter layer hopping becomes more like diffusive. Hence a 
crossover into a non-metallic behaviour occurs. 

Note that there is a crucial separation of scales. We 
assume that, because of the highly anisotropic elec- 
tronic dispersion in these quasi-2d materials, £^7™™° ^ 
Tboson 3> t c . Thus, it is possible for the strong electron- 
boson coupling to show up dramatically in c-axis trans- 
port, while the in-plane behavior is relatively unaffected. 
Also, it allows us to treat the c-axis hopping pertur- 
batively and still have quantitative control over the 
crossover region near Tb OS on- 

In this paper, we shall concentrate on the d.c. resistiv- 
ity in the c-direction. Our main finding is, indeed, that 
due to the small polaron, a metal to non-metal crossover 
for temperatures near the boson scale Tb OS on is obtained. 
At higher or lower temperatures, the c-axis transport 
basically reflects the in-plane dynamics through the in- 
plane scattering rate. 

We have so far not specified the nature of the bo- 
son. This is because the analysis to be presented can, in 
general, be applied to any kind of strong electron-boson 
coupling, provided there is this separation of scales men- 
tioned, and the bosonic mode is gapped and neutral. For 
the rest of the paper, we shall phrase the analysis in terms 
of phonons, with a view towards application to some 
members of the strontium ruthenate family. For other 
members of the same family with known magnetic insta- 
bilities, and perhaps for the cuprates, this boson could 
be a paramagnon or other magnetic excitation. 

Strontium ruthenate Sr2Ru04 is iso-structural to one 
of the parent compounds of the cuprates, La2CuC>4. 
Above its superconducting transition temperature of 
1.5K, it is a well-characterised 3d Fermi liquicO, and 
yet, above a coherence scale of ~25K, the c-axis optical 
conductivity loses its Drude peaklj, and above a "cross- 
over" scale, T max ~ 130K, the c-axis resistivity becomes 
non-metallic while the a6-plane resistivity stays metallic, 
albeit with a close to linear temperature dependence and 
a value which eventually exceeds the MotUloffe-Regel 
limha. In the recent report of Jin et. alS3 in the re- 
lated system Cai.7Sr .3RuO4, where there is the same 
dichotomous resistivity behaviour, T max ~ WOK appears 
to be tied to a structural phase transition. In Sr2Ru04, 
there is no such direct link, but there ate, strong electron- 
phonon interactions in the c-direction.119 This motivates 
us to consider phonon-assisted hopping as a potential ori- 
gin for the non-metallic p c in these systems at high tem- 
peratures. 

A further experimental observation in the vicinity of 
Tmax in the c-axis resistivity is that the magnetoresis- 
tance can change sign and becomes negative above T max . 
This happens^ in both the longitudinal (field parallel to 
the current) and transverse (field perpendicular to the 



current) directions in Sr2Ru04. Somewhat similar prop- 
erties can be found in some of the cuprates, although 
there seems to be no universality in the c-axis magnetoce* 
sistance between the various families of the cuprates.EEl 
In this paper, we shall calculate the orbital component of 
the magnetoresistance. This will necessarily be positive 
and will not explain these observations. However, a more 
complete theory would need to take into account the ef- 
fect of magnetic field on the propagators and the bosonic 
mode, which is beyond the scope of this paper, but we 
will argue that it must be these contributions that change 
the sign of the magnetoresistance, within this picture at 
least. 

The plan of the paper is as follows: in Section II, we 
define our model, in Section III, we study the model in 
the strong electron-phonon coupling limit. Section IV 
gives the results. Section V discusses the relevance of this 
model to the strontium ruthenate family and beyond. 

II. MODEL 

As already emphasized in the Introduction, the anal- 
ysis to be presented can be applied generally, with little 
modification, to any strong coupling of electrons to a 
gapped, neutral bosonic mode, but to be concrete, we 
shall consider coupling to phonons. 

Our model involves electrons coupled to the polar- 
isation induced by the out of phase vibration of the 
ions within a unit cell; the form of this electron-optical 
phouan coupling has been derived long ago in the 3d 
caseEil: 

^e-ph =Y1 I d 2 xY,Mq e M^Rn)Pn (x) (o ? +oL ? ) ,(3) 
n J q 

where R n = [x, z n ] T , x is the intraplane coordinate, 
z n = nc is the plane coordinate, with c the interlayer dis- 
tance and n the plane index. p n {x) — cj 1 (x)c„(a;) is the 
electronic density. The electrons mainly couple to lon- 
gitudinal phonons (represented by a^, at) because only 
these phonoiaa set up strong enough electric fields when 
they vibrateE3. Then the direction of the polarisation 
vector locally is simply along the direction given by the 
vector sum of the displacements of the ions in a unit cell. 
Now, for the mode that can couple to c-axis charge trans- 
port, we expect that in these highly anisotropic metals 
such as the cuprates and Sr2Ru04, the effective force 
constant between the planes may be less stiff than those 
within the plane. This is because of the strong cova- 
lent bonding and the strong band electronic correlation 
within the plane, as compared to the weaker bonding and 
correlations between the planes. Effectively, due to the 
large in-plane effective force constants, we let the vibra- 
tion frequency for the in-plane modes ujq* — > oo com- 
pared to the interplanar one. In any case, we wish to 
separate out the direct effect of strong electron-phonon 
coupling on c-axis transport, from the effect of phonons 
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on in-plane propagation which then influences the o 
axis transport via the mechanism sketched in Fig. [I]. 
Thus, with such a simplification, the Fourier transform 
of the displacement vector Q is only in the interplane 
direction: Q oc f (w ? ) _1 / 2 (a f + oL f ) w [0,0,Q C ] T , with 

Q c oc (L) q )~ 1 l 2 (a q + aL„), and the propagation direction 
q (not a vector!) is along the c-axis. In effect, all the ions 
in a plane oscillate together rigidly, in the direction per- 
pendicular to the plane, and the propagation direction is 
also parallel to the c-axis. Hence it is the total charge 
Pn = J d 2 xc^fflcnlx) of the n'th plane that couples to 
the phonon creation and annihilation operators a q ,a q : 



H c - P h = >J M q exp(iqz n )pr, 



and also we are able to obtain many explicit analytical 
results. We shall also discuss somewhat more qualita- 
tively the opposite case of a generic dispersion where the 
phonon density of state does not have any sharp features. 
These two forms of dispersion can be thought of as op- 
posite limits: the Einstein phonon is the limit where all 
phonon modes (of different wavevector q) have the same 
energy, while the more general dispersion corresponds to 
the limit when all phonon modes are non-degenerate. 



III. METHODS 



(4) 



Canonical Transformation 



where we have combined all the prefactors as usual into 
the clcctron-phonon coupling M q . 

If the bosonic mode is of magnetic origin, we may ar- 
gue that because of stronger electronic overlaps in-plane 
compared to the interplane direction, and hence stronger 
in-plane exchange effects, the magnetic fluctuation mode 
is softer in the interplane direction, and again this sim- 
plifed form of electron-boson coupling will be appropri- 
ate. 

Our model Hamiltonian is thus: 



H 



J2H {n) +H c + H e _ ph + H ph , 

n 

/ d 2 x^^t c c s n+1 (x)c n (x) + H.c. , 
J n 

M q exp(iqz n )p n (a q + al g ) , 

n q 



The model system consists of a stack of 2d planes de- 
scribed by H^ n > , We shall take this as a phenomeno- 
logical input, characterised by a 2d Green's function or 
equivalently, a 2d spectral weight. 

In H c , electrons are assumed to hop from one plane 
directly to a neighbouring plane only, at the same x; in 
momentum space, this means that the c-axis hopping ma- 
trix element t c is independent of the in-plane momentum. 
(This ignores the complication of staggered planes in the 
cubic perovskite structure and the multiband nature of 
Sr2Ru04, see Discussion.) For simplicity, we have left 
out the spin index on the electron operators c n {x). 

-He-ph describes the electron-optical phonon interac- 
tion in the c-direction, as detailed above. We have left the 
electron-phonon coupling M q unspecified: as we shall see, 
it enters only within an electron-phonon parameter [A or 
7(T) depending on the phonon dispersion, see Methods]. 

For -ffph, The c-axis optical phonons have a dispersion 
Lj q \ in this paper we shall look closely at the Einstein 
phonon oj q ~ luq. This is both because optical phonons 
tend to have little dispersion relative to the acoutic ones, 



The key physics we are looking for is the effect of strong 
electron-phonon coupling on the charge transport in the 
weakest direction, the c-axis. Physically, the motion of 
the electron is accompanied by the emission and absorp- 
tion of a large number of phonons due to Xkc strong cou- 
pling, forming the so-called small polaront3. Technically, 
one must deal with the large term 7J c _ p h first, and add 
the inter-plane tunneling H c later as a perturbation. The 
former is accomplished by the canonical transformation, 
H = exp(—S)Hexp(S), a straightforward generalisation 
of the transformation well-known in the small polaron 
problemE3'Ej: 



S 



d xS(x), 

l0„ 



S(x) = -^2^e^cUx)c n (x)(a q -al q ). (6) 
One finds then, 



it . q 



H 



HAx) 



X„ — 



J2H^+ f d 2 xH c (x) + H ph , 

J2t c cl +1 (x)c n {x)xl +1 X n + H.c. , 

n 

-p{E^ (», -'-,)}. 



(7) 



Thus strong electron-phonon coupling leads to, (I) a 
renormalisation of the in-plane chemical potential (first 
term of H c ), which we shall henceforth ignore, and (2) 
an effective vertex correction for the c-axis hopping t c 
(second term of H c ). Note that since the phonon modes 
couple to the total charge density in a plane, H^ n ' — 
exp(—S)H( n,) exp(S') = H^ n >: the in-plane motion of the 
electron is not in the presence of the interplane polaronic 
effects, within this simplifed model. 
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B. Linear Response: conductivity and 
magneto-conductivity 

In this paper, we look at the zero- frequency con- 
ductivity and the magnetoconductivity. For the con- 
ductivity, since the charge in the n'th plane is Gf n = 
e J d 2 xc} l (x)c n (x), the current in the c-direction is just 
3n — dtQ n = —i[Q n , H]. After the canonical transforma- 
tion, j n — ~i[Q c nl H c ], and so: 



tin) = i^c 



d 2 x 



C n+l(x) C n(x)Xl +l X. 



(8) 



Using the Kubo formula and expanding to 0(t 2 ) gives 
the conductivity, ie., the linear repsonse to an applied 
electric field in the c-direction a c = j c /E c : 

e 2 t 2 f f 
a c (k c = 0, ft = 0) = lim ^ / d 2 k\\ / dr (9) 

n J J 

U(r) = (xt +1 (r)X n (T)Xt(0)X n+1 (0)) . 

where Gn d \k\\) is the in-plane (dressed) Green's function 
for the plane n at in-plane momentum /c|| . In Fourier 
space, the effect of strong electron-phonon coupling is to 
induce an effective w-dependent vertex t 2 — > t 2 U(w). 
Classically, the phonon correlation function U[ui) leads 
to a Debye- Waller factor. 

For magnetoresistance, we shall only study the case 
where the magnetic field lies in-plane, avoiding the com- 
plication of magnetic field induced orbital effects on the 
in-plane propagators. We also do not consider spin ef- 
fects here, see Discussion. As B does not couple directly 
to the phonons, its sole effect is on the orbital motion 
of the electrons, which can be dealt with by the usual 
Peierls substitution p — ► p — eA. With B = (5,0,0), 
the associated vector potential (in the Landau gauge) is 
A = (0, — zB, 0). Choose plane n = 1 to have c-axis co- 
ordinate z = 0, and for plane n = 2, z — c where c is the 
interplane distance. Hence, 
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Gn = l(fc||), 

ecx B, 



(10) 
(11) 



where c is a vector of length c pointing in the positive 
c-direction, and fc|| is the in-plane momentum. 

To massage the (magneto) conductivity into a more 
familiar form, we introduce the usual spectral represen- 
tations: 



G^Hk 



dz 



A^ 2d \k h z) 



U{v n ) = I dz- 



lUJ m - z 

B{z) 



where uo m = (2m + \ )nT is an odd Matsubara frequency, 
and v n — 2nirT is an even Matsubara frequency. T is 
the temperature. Note that since there is translation 
invariance in the c-direction, we have dropped the plane 
index n in the spectral functions A^ 2d ^(k^ 7 z) and B(z). 
Fourier transforming with respect to the imaginary time 
in Eq. (^), and doing the Matsubara sums leads to the 
d.c. conductivity: 



e 2 N c ti 



dio dv 



d 2 h 



T J J (2tt) 2 

A<M{k\\ +q B ,uj + v)A ( > 2d \k h v) 
j{u)-f(Lo + u) 



(12) 



D(ui,T) = lu B(u)n"(u) [I + 



where N c is the number of planes in the c-direction, 
f(v) — l/(exp[z//T] + 1) is the Fermi function, and 
n B (uo) = l/(exp[w/T] - 1) is the Bose function. D(u,T) 
is essentially the phonon spectral weight multiplied by 
the Bose factors. 



C. Phonon Spectral Function 



To get the phonon spectral function B(to), we need: 



B(uj) = Im(7 ret H 

7T 



-— Im 

7T 



e Mt [/ ret (i),(13) 



where the retarded real time spectral function is just: 
U rct (t) = -i6(t)(\xl +1 (t)X n (t) , At(0)A„ +1 (0)j ) (14) 

with X n (t) given in Eq. (?]. 

The calculation of U let (t) proceeds along Xhe same line 
as in the small polaron model (see MaharO, Ch. 6.2), 
and we get the exact result: 

U Ict (t) = -i9(t)2e^i Fq[nB WUW 5 ^K))-( 2 « B 'K)+i)l 

x {cos [oj q {t + i/2T)} - cos [u q (-t + i/2T)]} (15) 



2(1 — cosq) > 0, and n B (tu q ) 



where F q 

l/(exp( Wg /T) - 1): 

In this paper, we look at two types of phonon disper- 
sion: (1) Einstein phonon where tu q = loq for all q, and (2) 
a generic dispersion such that there are no sharp features 
in the phonon density of state. Again, the calculation is 
directly analogous to tha-Lof the small polaron problem, 
see e.g. chapter 4 of Refill for details of the calculation. 

For the Einstein phonon, setting co q = ljq in Eq. [l^, 
putting this form into Eq. O, and using the Bessel func- 
tion identity exp(z cos((/>)) ^^Z^Tloo ^n(z) exp(in<j)), the 
spectral function can be expressed as an infinite sum of 
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Bessel functions of the imaginary kind: 

B(u) = 2e~ 2ST (16) 

x £ J n (A(T))sinh(^) 5(0,-^0), 



Sr 2 Ru0 4 , a = 3.87A c = 6.37A around room tempera- 
ture, and so, 



to 



3-10" 4 -B (5intesla). 



(22) 



n— — 00 



where 

A(T) = (A/ Wo ) 2 /sinh(^ /2T), (17) 
A 2 = ^|A/ 9 | 2 2(l-cos ), 

25 T = (A/w ) 2 coth(w /2T). 
Thus the function -D(w, T) defined in Eq. O becomes: 



D(u,T)=e- 2ST ^TI (A)5(u) 



Thus for the in-plane electronic dispersion, we can ap- 
proximate — (-k + V • <fs with = de^/dk. 

We simplify further with the choice of a flat band of 
width 2D and a circular (in-plane) Fermi surface. Then 

/ d 2 k — > No J° D defjf* |§, where N is the 2d density 
of state at the Fermi surface. For the isotropic in-plane 
Fermi surface, we pick q B to lie in the direction of k Xl and 
thus, e^r + - s — fife + Uf<7b cos(0). For simplicity, we have 

taken B to lie entirely in the plane, hence q B = ecB. 
Write the 2d spectral weight in the form: 



(18) 



A^(k, 



2J (A) [5(uj - nw ) + 5(w + nw )] 



1 



7T (<J 



r(«) 



(23) 



2sinh(w/2T) ^ 

For the more general dispersion, the time integral has 
to be approximated. As in the small polaron problem, 
one appeals to the fact that this is at strong electron- 
phonon coupling, and so the dimensionless factor F q in 
the exponential in Eq. ^| is large, and the integral can 
be done using the saddle point approximation: 



B(u) 



7 ^sinh(w/2T)e- w2/472 



(19) 



ex p \ - F i \f nB ( uj i) + i - \A^K) 



with 7 2 = Y^ q \ M q\ 2 ( l - cos g)/ sinh(w g /2T) and F q is 
as before (see after Eq. |l^). And the function D(to,T) 
becomes: 



where we assume that the scattering rate T(u)) has little 
in-plane momentum dependence. As all energy scales 
will be small compared to the bandwidth 2D, we can set 
D > 00, and we find then, 



d 2 kA (2d) (k + q B ,Lo + v)A {2d) (k, to) (24) 
^Re [(v F q B ) 2 - {u + iTtotCT, w, v)) 2 } - 1 ' 2 , 



where T tot (T,u,v) = T(T,lj) + T{T,u; + v). 

Substituting Eq. |24| into Eq. [l2] leads finally to the 
zero-frequency, zero-momentum c-axis conductivity: 

<7 C (T, B) = + 



ttT 



D{u,T) 



21V 



D(u,T) 



1 



2 /47 2 



4y^2 sinh(w/2T) 



D. Calculation of a c 



[(v F q B ) 2 - (f + iT to t(T, to, v)) 



This is the key equation we shall analyse in the next 
section. It is worth summarizing the assumptions used 
in deriving this: (i) B lies in-plane and is small, (ii) flat 
band, circular Fermi surface, (iii) in-plane spectral weight 
has a scattering rate that has little momentum depen- 
dence, and (iv) second order perturbation theory in t c . 



It is advantageous to do the in-plane momentum inte- 
grals first. Note that the ratio of the magnetic wavector 
q B = ecx B (Eq. |Tl| ) to the in-plane Fermi momentum 
k F is: 

qs ^ ecB _ 2acB 

k F ~ n/2a ~ $ ^ ' 

where $0 = h/2e is the unit flux quantum, a is the in- 
plane lattice spacing, c is the interlayer distance. For 



IV. RESULTS 

A. Einstein phonons 

Because of the totally degenerate spectrum uj q — luq, 
the function D(tu,T) (Eq. |l8|) is made up of delta funtions 
at the harmonics u) = najQ, n = 0, ±1, ±2, . . .. Putting 
this form into Eq. ^5] for a general in-plane scattering 
r(T, uj), we get the exact result: 
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e 2 N N c t 2 C c _ 2ST 



TV 
OO r 



-df 



oo(A(T)) 



n=l 



a„(A(T)) Rc 



^ v/(« F e C B) 2 + (r tot (T,0^)) 2 

1 



(26) 



a/ (v F ecB) 2 + (T tot {T, nuj , v) - inuj )' 



+ {ujq — * -w ) 



where a (A) = I (A),a n (A) = In(A) Bin ^°/,j 2T) , and 
I n {z) are the Bessel functions of imaginary kind. One can 
show that the amplitude, et„, decreases with increasing 
n. In the experimentally accessible range of temperatures 
(at most a few times of Tuebye ~ 500K), one only needs 
the first few harmonics (n) in Eq. (|26]). This is because 
the higher harmonics contribute little to the low T and 
small B behaviour, thanks to both the amplitude a n (A) 
getting smaller and the denominator of Eq. getting 
larger with higher n. 

To determine the asymptotes, note that the various pa- 
rameters of the Einstein phonon model have the limiting 
values (see after Eq. Fyl for their definitions): 



A(T) 



(wo) 



A) 2e~ UJo/2T 



T«w , 



2S* T 



2T 

UJQ 

2 



T » w , 
T « w , 



(27) 



where A 2 = J2 q |A^| 2 2(1 — cosg) characterises the 
strength of the electron-phonon interaction. The argu- 
ment A(T) of the Bessel function is always small in the 
low T limit, but it can be large or small in the high T 
limit, depending on the magnitude of A/wo via the com- 
bination (A/w ) 2 (T/w ). 

Substituting the asymptotic forms Eq. ^ into Eq. |2^, 
and using the expansions for Bessel function, we get 
for low temperatures T <C ujq, (or for high T with 
(A/w ) 2 (7> ) « 1): 



e 2 N N c t 2 c 



-(A/wo) 2 



^(v F ecB) 2 + 4r(T,0) 5 



(28) 



We have dropped n > terms, as they are of order 
e —nuj /T re i a tive to the n = term. Also, both the 
marginal Fermi liquid and the Fermi liquid scattering 
rate varies relatively slowly compared to df/dv, and 
hence we have approximated —df/dv ~ 8{v) to get the 
above equation. Note that in this low T limit, the effec- 
tive bandwidth t c gets renormalised by the phonons by 
the factor exp[-(A/w ) 2 /2]. 



For high temperatures T ^> ujq and A/luq > 1: 
e 2 N N c t 2 c l + 0{u Q /T) 



7T 




< A S 


2 




\ 2T 








/ "0 


1 


1 







(29) 



y/(v F ecB) 2 + 4r(T, 0) 2 
1 



+2Rc 



x /(vFecB) 2 + (2r(T, W o)- J w ) 2 



+ (higher harmonics) j . 

Again, we have made the approximation (/(f) — f(v + 
nu>o))/nLOo ~ —df/dv ~ <5(f), the first approximation 
being valid as T > wo, and the second approximation is 
due to the same reason as mentioned for the low tem- 
perature case. Notice that the leading terms in both the 
high- and low-temperature limits do not have any e~ A//T 
factors, in contrast to the general phonon dispersion case 
below. 

The precise low and high T- dependence at B = re- 
flects the form of the in-plane scattering. For a marginal 
Fermi liquid form: 



T(T, uS) — (»mflw coth (u>/T) . 

the asymptotics from Eqns.|2^,|2^ are: 

ct c (T<w ) ~ T- 1 
a c (T>w ) ~ T~ 3/2 

while for a Fermi liquid form: 

u? + (tt T) 2 



(30) 
(31) 



r(T,u>) = a FL - 

the asymptotics are: 

cr c (T < w ) 
cr c (T > w ) 



Tl — 2 

ji-5/2 



(32) 



(33) 



The above asymptotics miss out the experimentally 
important and interesting regime of T «s Wo- A broad 
maximum appears when T is somewhat below u>o, only 
when A/ujq > 2. Also, the maximum grows rapidly with 
increasing A/luq. This occurs because of two compet- 
ing factors: see Eq. p6l For A/loq 3> 1, A(T) is again 
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large, and we can use the asymptotic form for the Bessel 
functions, while keeping T « lu : 



sinh 2 (o; /4T)\ 



a c oc y sinh(wo/2T) exp 
f 1 



- 


'A" 




v _ 



h(^ /2T) ; 



[(7; F ecS) 2 + (2r(T, 0)) 2 ] 1/2 



(34) 



Fig. H illustrates the resistivity (the inverse of the con- 
ductivity) of Eq. ^6] for a range of values of A, using a 
marginal Fermi liquid form for the in-plane scattering, 
and including up to n = 5 harmonics of wo- Fig. [| shows 
the Fermi liquid case, where the maximum is much less 
prominent than ip_the marginal Fermi liquid case, as al- 
ready mentioned.EJ All energies are measured in units of 
loq. p c is measured in units of ir/e 2 NoN c tl. 



where we have dropped higher harmonics (which change 
the magnitude, but not the general T-dependence) . Thus 
there are two opposing tendencies: the phonon con- 
tributed part (the exponential and the square root of 
sinh) which grows with T, versus the in-plane scatter- 
ing part r(T, 0) oc T for the marginal Fermi liquid. The 
result is the broad maximum structure in er c (Fig. 
Note that the height of the maximum has a strong de- 
pendence on the strength of the electron-phonon coupling 
A. In a Fermi liquid the in-plane scattering rate ~ T 2 
damps out the phonon contributed parts more strongly 
and to see this structure, one needs to go to larger A/wo- 
We have estimated the critical phonon parameter 
A c /u)q and the position of the maximum T max and the 
concomitant minimum T m j n from Eq. ^6] by taking into 
account only the zeroth harmonics n — 0. For an in-plane 
scattering rate of the form T(T) = aT v , 



(35) 
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FIG. 2: B — c-axis resistivity with in-plane marginal Fermi 
liquid scattering (Eq. ^) for electrons coupled to c-axis Ein- 
stein phonons, qmfl = 1.0 



and the broad maximum only appears when A > A c . 
The critical temperature when the broad maximum 
starts to appear, i.e. when T max starts to be different 
from T m i n , occurs at: 



T, 



crit 



■■ wq/3.0. 



(36) 



For A only slightly larger than A c , both T max and T n 
are close to T cr it: 



w 



2T, 



max 
min 



^0 



± 



\ 



V 

0.2 



^0 \ _ /woV 
A J V A J 



(37) 



For large A/oj , T n 
equations: 



and T m j n obey the approximate 



A / C0 = 2.2 




T/ co. 







A 2 



A 2 



I1K1N / W 

exp 



V2T n . 



(38) 



L0 



Note that both T max and A c depend mainly on phonon 
parameters; one can show that the in-plane scattering 
rate enters only in the form of the exponent rj. In partic- 
ular, both T max and T m ; n can only depend on the scale 
luq. Hence the width of the broad maximum, roughly 



given by 2(T lr 



is governed by the scale loq. 



FIG. 3: B — c-axis resistivity with in-plane Fermi liquid 
scattering (Eq. |32| ) for electrons coupled to c-axis Einstein 
phonons, ofl = 5.0, £f = lOOt^o. 

Now we consider the magnetic field dependence. The 
magnetoresistivity is defined as Ap c — p c (B)/p c (0) — 1. 
For later convenience, the effective cyclotron frequency 
is defined as: io c — v^ecB . 

First, at very low T and very large magnetic fields, 
in Eq. the effective cyclotron frequency can become 
larger than the scattering rate and so the magnetoresis- 
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FIG. 4: c-axis resistivity at finite B with in-plane marginal 
Fermi liquid scattering (Eq. ^) , for electrons coupled to c-axis 
Einstein phonons as a function of T. «mfl = 1.0, A/uo = 2.0. 
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FIG. 5: c-axis resistivity at finite B with in-plane Fermi liquid 
scattering (Eq. for electrons coupled to c-axis Einstein 
phonons, as a function of T. ofl = 5.0, £f = lOOtdo, A/uio = 
2.0. 



While it is possible for the coefficient of the B 2 term in 
Eq. |o| to chan ge sign, it can be shown that the n = 
term always dominate, leading to a 1 — B 2 conductivity, 
or a +B 2 magnetoresistance: 

A Pc cx (w c /T cff ) 2 , (41) 

where the effective scattering rate r e ff is made up of both 
the electronic scattering and the phonon frequency, and 
taking into account only up to n — 1 harmonics, 

1 , 4T_i 

cff 1 , 2r 1 (4r;-3o J g) ' K ' 

Terf + (4Tf+^? 

where T = T(T,0), and Ti = T(T,u ). 

In conclusion, there is positive magnetoresistance at 
all temperatures and for both Fermi liquid and marginal 
Fermi liquid in-plane states. In general, if the in-plane 
scattering T(uj) could be made larger at u> — relative to 
the higher harmonics u> — ntoo, the contribution of the 
n = part of the phonon function D(u>) is suppressed, 
and this can lead to a negative magnetoresistance. 

Figs. and ^ illustrate the resistivity as a function 
of T and B. 



60 
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tance goes as: 



Ap c oc w c /T(T), 



(39) 



i.e., linearl$B in B. See Fig. 

At higher T, the scattering rate dominates in both 



Eqs.||and34. For the n = term, [w 2 +r tot 2 ] -1 / 2 gives 



a 1— B 2 contribution to the conductivity, i.e. a +B 2 mag- 
netoresistance. For the n > 1 terms, since luq ~ 10 _2 ep, 
and w c ~ 10 _3 eF for B up to a few tesla, u> c -C Wo, wc 
can extract the leading B behaviour of the n > 1 terms 
in the conductivity: 



Re- 



1 



a/w 2 + (r to t - mw ) 2 

3(n^o) 2 



F 2 

1 tot 



r 2 

1 tot 



2 [r 2 ot + (™ ) 2 



+ (noj ) 2 
+ O(cj c /T tot ) 4 



(40) 



FIG. 6: c-axis resistivity with in-plane marginal Fermi liquid 
scattering (Eq. for electrons coupled to c-axis Einstein 
phonons, as a function of B. omfl = 1.0, A/luq — 2.0. T is 
measured in units of ujq. 



B. General phonon dispersion 

Having considered an Einstein mode, we now treat 
a more general dispersing phonon mode. The func- 
tion D(uj,T) will be seen to be strongly peaked near 
lu ~ 0, compared to the in-plane Fermi energy scale cf- 
Furthermore, we are in the regime T/eF *C 1. Thus, 
[f(v) — f{v + uj)]/uj ~ —df/dv ~ 5(i^). The last approx- 
imation is valid for the same reason as in the Einstein 
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FIG. 7: c-axis resistivity with in-plane Fermi liquid scattering 
(Eq. |3^), for electrons coupled to c-axis Einstein phonons, as 
a function of B. ofl = 5.0, £f = lOOc^o, A/ljo = 2.0. T is 
measured in units of loq. 



phonon case. Hence, with these approximations and in- 
serting the gaussian form of the phonon function of Eq. |2^ 
into Eq. GB leads to 



e 2 N a N c t 2 c 
irT 

xRe 



-C(T) 



dv 



C(T) 



V^(T) smh # 

ex P ~4 7 a (T) 

[{v F ecBf -{u + iY tot {T,v)f] 1/2 



\ AT I 



sinh ^ 



where F„ 



M„ 



2(1 - cosg) and 7 2 (T) = 

Y Jq F q \u q \ 2 /2s\nh{Lo q /2T). 

The various parameters have the following limiting val- 
ues: 



7 2 (T) 
cxp -C(T) 



2TA : T » ui q , 
7o 2 (T) : T« Wg , 

g— A/T . T>w 



?! 



(44) 



A = J2 q F q Lu q /2 involves the electron-phonon coupling, 
together with the phonon dispersion. It characterises 
the strength of the electron-phonon interaction, and is 
directly analogous to the parameter A for the Einstein 
phonon case. 7o(T) = Y2q F q \to q \ 2 exp -u q /2T. (Note 
that the expressions T 3> uj q and T <C u q refer to some 
typical energy of the phonon dispersion, e.g. the Debye 
frequency.) 

In the low temperature limit T -C Lu q , the phonon 
spectral function provides a gaussian in frequency with 
a width 7o(T) which is exponentially narrow in T. Thus 
in the integral over frequency ( J dv) , only v rj con- 
tributes. Furthermore, since T 3> 70 (T), the scattering 



rate T(T, v) can be taken at zero frequency, for both a 
Fermi liquid and a marginal Fermi liquid. Then the in- 
tegral in Eq. [l3] can be done immediately, leading to the 
low T asymptotic form: 



a c (T< Uq ,B) 



2e 2 N N c r c e~^* 
TTy/{v F ecB) 2 + 4r(T,0) 2 



(45) 



This gives rise then to a low temperature metallic 
resistivity that tracks the in-plane scattering T(T, 0), 
and a positive magnetoresistance with a leading in- 
dependence, much as in the Einstein phonon case. Note 
that the hopping matrix element t c is effectively reduced 
by the phonon-induced prefactor exp — ^ q F q /2, similar 
to the Einstein phonon case. 

In the high temperature limit T 3> uj q , when A <C T, 
the width of the gaussian 7 2 ~ 2TA <C T 2 , hence again 
only small frequency v ~ is needed in the integral in 
Eq. W3, leading to: 



a c (T^uj q ,B) 



2e 2 N Q N c t 2 c e 



-A/T 



TTy/(v F ecB) 2 + 4r(T,0) 5 



(46) 




A / CO = 2.0 
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FIG. 8: B=0 c-axis resistivity with in-plane marginal Fermi 
liquid scattering , for the phonon dispersion uiq = uja\l — 
0.4cos(q)], qmfl = 1.0. 

Just as in the Einstein phonon case, there can be 
a broad maximum structure when T ~ uo q : it can 
be shown that in this regime, the phonon contribution 
e -c(T) /^(j 1 ) provides a semi-conducting-like upturn in 
conductivity, which competes with the in-plane scat- 
tering r(T, 0) that eventually dominates as T » uj q . 
This is to some extent captured in Eq. above. More 
precise T-dependence will require an explicit form for 
the optical phonon dispersion tu q . For illustration, let 
Lo q = uj {1 — 0.4cos(g)], and set the electron-phonon cou- 
pling M q to be ^-independent. Eq. ^ is then evalu- 
ated numerically and plotted: Fig. |s| is for a marginal 
Fermi liquid form for the in-plane scattering, while Fig. ^ 
shows the Fermi liquid case. p c is measured in units of 
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FIG. 9: B=0 c-axis resistivity with in-plane Fermi liquid scat- 
tering , for the phonon dispersion u q = wo[l — 0.4cos(g)], 
«fl = 5.0, ef = lOOwo • 

Tr/e 2 NoN c tc- Qualtitatively, these plots are similar to 
the Einstein phonon ones. 

As for the Einstein phonon case, there is only positive 
magnetoresistance for the more general phonon disper- 
sion, at all temperatures. 

V. DISCUSSION 
A. Summary of results 

In this paper we have studied 2d planes of electrons 
characterised by a Fermi liquid or a marginal Fermi liquid 
scattering rate (that has little momentum dependence), 
coupled by a small inter-plane single-particle hopping in 
the third dimension, in the presence of a strong inter- 
action between electronic density and a bosonic mode 
in the inter-plane direction. The electron-boson cou- 
pling is treated exactly via a canonical transformation, 
and we treat the inter-plane hopping to lowest order (i 2 .) 
to calculate conductivity and magnetoconductivity. For 
concreteness, we have considered the bosonic mode to 
be a phonon, and have treated two types of phonons: 
Einstein phonons and a more general dispersing phonon 
mode. The overall behavior appears to be quite similar 
for the Einstein mode when compared to the more general 
phonon dispersion, and we shall discuss them together, 
pointing out the differences where they arise. 

We find that at temperatures low compared to the typ- 
ical phonon energy, the c-axis resistivity p c is metallic, 
displaying the same temperature dependence as that of 
the underlying in-plane scattering rate T(T): p c ~ T(T). 
The magnitude is however strongly reduced, correspond- 
ing to a much smaller effective c-axis bandwidth, directly 
analogous to the bandwidth rjeuormalisation in the Hol- 
stein small polaron problemE£lEj. 

At temperatures much much larger than the typical 
phonon energy, (a regime that is probably academic,) 



p c (T) again reflects the T-dependence of the in-plane 
scattering rate, but in the Einstein phonon case, there 
is an extra multiplicative factor of T 1 / 2 contributed by 
the phonons, while for the general phonon case, the T- 
dependence is exactly as in F(T). Note that it is because 
of this in-plane factor of T(T, u>), that the high tempera- 
ture state is metallic like {i.e. dp c /dT > 0), in contrast 
to the standard Holstein small polaron model. 

In the experimentally interesting regime where T is of 
the same magnitude as the phonon energy (a few hundred 
Kelvins), we find the possibility in p c of a broad maxi- 
mum structure that connects the low and high tempera- 
ture metallic states. For the Einstein phonons, we have 
shown that this broad maximum, i.e. a "non-metallic" 
T-dependence, occurs only when there is strong electron- 
phonon coupling (Eq. [35]). We have also estimated the 
position and width of the maximum in the temperature 
axis. For the general phonon dispersion case, we expect, 
and see numerically, qualitatively similar behavior. Note 
that the physics is richer here than in the Holstein small 
polaron problem, due to the extra degree of freedom of 
the in-plane scattering rate. 

As expected, we have found only positive magnetore- 
sistance in our model, since we only consider the effect 
of the magnetic field pointing in the plane on the orbital 
motion of the electrons. This goes as Ap c oc (w c /T e ff) 2 
(Eq. |4l]), except for the regime when the magnetic field 
is vexy large and the temperature is very low, in which 
caseEa Ap c oc oj c /T(T) (Eq. ^). This is so for both the 
Einstein phonons and the more general phonon disper- 
sion, and for both the marginal Fermi liquid or Fermi 
liquid scattering rates. 

In this work we have tacitly assumed that comput- 
ing o~ c as an expansion in t c is valid and that we can 
fully characterize the in-plane physics by its single par- 
ticle properties. For non-Fermi liquid states in the plane 
this may not be.-fclie case (a Luttinger liquid in ID is a 
notable exampleEa). Under such circumstances our re- 
sults would presumably breakdown below a temperature 
scale which would depend on the details of the in-plane 
physics. 

Our model consists of phonons that have no struc- 
ture in the in-plane direction, and a more realistic model 
should have a fully three dimensional dispersion. In this 
respect, a question immediately arises: if the metallic 
to non-metallic transition in the c direction requires a 
strong electron-phonon coupling, might not this strong 
coupling also show up in the in-plane physics? This 
is presumably a question of scales: the characteristic 
phonon energy is expected to be larger than the c-axis 
hopping matrix element t c , but the phonon energy is still 
much smaller than the in-plane Fermi energy, because 
of the highly anisotropic nature of the materials consid- 
ered here. Hence, one expects the strong electron-phonon 
coupling to have a much more dramatic effect in the c- 
direction, than in the in-plane direction. 
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B. Application to the ruthenate systems 

There exist optical phonqas with the appropriate sym- 
metry for c-axis transporter in S^RuCU, and experi- 
mentally the broad maximum in the c-axis resistivity 
has been linked to a structural phase transitionlij in 
Cai.7Sr .3RuO4 at around the broad maximum tempera- 
ture. Thus one should consider the possibility of electron- 
phonon interaction affecting the c-axis transport. 

In particularjJbpth the Sr 2 Ru04 and Cai.7Sr .3RuO4 
systems exhibittfllij qualitatively this broad maximum 
structure found in our simple model, in the c-axis resis- 
tivity near to their characteristic (c-axis) phonon energy. 
In Sr2RuC>4 there is even tha-iint of an upturn in re- 
sistivity at high temperaturesEH, and the result plotted 
in Fig, H bears some resemblance to the data of Tyler 
et.a/liil (their Fig. 3). Also, the broad maximum temper- 
ature of ~ 130K in Sr2RuC>4 is indeed smaller than [the 
characteristic phonon energy (around 500 to 800 K,ll3), 
consistent with our estimated bound T max < 3t^o- Using 
the value of 500K for uj , we estimate that the electron- 
phonon coupling parameter A is about two times of u)q. 

However, the negative magnetoresistance seen around 
the same temperature regime as the broad maximum can- 
not at present be explained within this phonon-assisted 
hopping picture. In our simple model here, to obtain 
negative magnetoresistance robustly, one needs either a 
phonon function D(u>) that peaks at a finite frequency, 
and/or an in-plane spectral weight where the scattering 
rate is smallest at a finite frequency. Either of these 
possibilities suppress the u> = contribution to er c (see 
Eq. |25|.) The Kondo effect can give rise to the second 
possibility, while it is unclear what kind of phonons can 
lead to the first possibility. 

Other factors that we have neglected that could in- 
fluence the magnetoresistance include the details of the 
Fermi surface. We have used a flat band with a circular 
Fermi surface, while experimentally, the (3 band (which 
has the largest c-axis dispersionliZI) is nearly a square 
Fermi surface and moreover, the c-axis electronic disper- 
sion has dependence on the direction in the in-plane mo- 
mentum space. Also, the Ru0 2 planes are staggered from 
one plane to another, so that an electron hopping from 
one plane to another does not stay at the same in-plane 
coordinate x as we have assumed in Eq. [s]. These facts 
may lead to a more singular momentum-integrated spec- 
tral weight, compared to the form in Eq. ^ used here. 
But if the spectral weight is now concentrated into a nar- 
rower range of momenta, the B-induced relative shift in 
the momenta between the planes would lead to a more 
positive magnetoresistance. 

Alternatively, there may be additional hopping chan- 
nels that contribute to the negative magnetoresistance, 
for example, hopping via some intermediate localized 
stateM 

Next, we discuss spin effects, which have been ignored 
so far. At least in the low temperature Fermi liquid 
regime (T < 25K), there are various types of spin fluc- 



tuations in JSr2Ru04 (both antiferromagneticEa and q- 
independcnt£3) . In-plane spin fluctuations enhance the 
in-plane electronic scattering rate. If the applied mag- 
netic field freezes these fluctations, then one can get 
a negative magnetoresistance, with the magnetic field 
pointing in any direction. However, this negative magne- 
toresistance will be seen in both p c and p a \y. In S^RuO^, 
c-axis magnetoresistance goes negative above T max for 
magnetic fields pointing in-plane or in the c-directiaa, 
but the in-plane magnetoresistance remains positivetj. 
Thus one can rule out suppression of in-plane scatter- 
ing due to magnetic field freezing out spin fluctuations. 
And if there were spin-fluctuation mediated c-axis hop- 
ping, magnetic field may again suppress this, leading to 
a positive magnetoresistance. 

After this work was completed, we learnt that the in- 
plane spectral weight at temperatures near and above the 
resisitivity maximum temperature may be anomalousEJ, 
in contrast to the spectral weight observed in the low 
temperature Fermi liquid state. In the Introduction, 
we have pointed to the possibility that anomalous in- 
plane physics can lead to incoherent c-axis transport, 
and these recent results suggest that for Sr 2 Ru04, this 
needs to be taken into account. However, because of 
the linkage of the structural phase transition and the re- 
sistivity maximum, the strong electron-phonon coupling 
model is still relevant to the anomalous c-axis transport 
in Cai.7Sro.3RuC>4. 



C. Boson-electron coupling 

Any strong boson-electron coupling can in general be 
tackled by the present technique. The analysis presented 
is for the specific case of a neutral, gapped, bosonic mode, 
as in a phonon. Instead of phonons, there could be 
magnons or paramagnons describing the effects of spin 
fluctuations. In addition to the cuprates, other mem- 
bers of the family of strontium/calcium ruthenate exhibit 
some types of spin fluctuations and work is in progress to 
check their importance in the c-axis transport. For mag- 
netoresistance the calculation however needs to be mod- 
ified. In this paper, the effect of magnetic field is purely 
on the orbital motion of in-plane electrons because there 
is no direct coupling between the magnetic field and the 
phonons, which is not the case for (para) magnons. 

D. Conclusion 

We have demonstrated that a strong coupling between 
electrons and a bosonic mode in the (weakest) c-direction 
in a highly anisotropic metal can lead to a broad max- 
imum (metallic to non-metallic crossover) in the c-axis 
resistivity. For the case where the boson is a phonon, we 
have discussed the potential application of the model to 
Sr2RuC>4 and its relative Cai.7Sro.3Ru04. Despite cer- 
tain simplifying features of the model, the qualitative 
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properties of this crossover in these layered ruthenates 
are captured succinctly. Work in progress will address 
the issue of the anomalous c-axis magnetoresistance in 
Sr2Ru04 within this framework, and also the effect of 
the anomalous in-plane spectral weight on c-axis trans- 
port in Sr2RuC>4. Also we are currently studying other 
signatures in transport properties of this strong electron- 
boson coupling, for example, in thermal conductivity. 

After completion of this work, we became aware of 
some related work by Lundin and McKenzie concerning 
a related small polaron model for intra- and interlayer 



transport [cond-mat/0211yyy]. We thank them for send- 
ing us a copy of their preprint prior to submission. 
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